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[57] ABSTRACT 
An improved brain wave electroencephalograph 
(EEG) method measures the physical positions of EEG 
electrodes applied to the subject's scalp. A set of stan- 
dard cross-directional lines are obtained to provide 
measurements of the subject's head and a computer 
system automatically defines the subject's head shape 
according to predetermined head shape classes. The 
position of the electrode on the subject's scalp is then 
determined by scaling the measurements of the subject's 1 
head* to a head model of the same head shape class. 

4 Claims, 5 Drawing Sheets 
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is measured. A flexible hat is placed on the subject s 
EEG SPATIAL PLACEMENT AND head and piezoelectric ribbons, embedded in the hat are 

ENHANCEMENT METHOD stretched. The ribbons are directed in cross-lines and 

the extent the ribbons are stretched (tensioned) deter- 
BACKGROUND OF THE INVENTION 5 mines the size of the subject's head in the direction 

1 F* M f I t* along which the ribbon lies. The head measurements are 
TheTnve^tionrela^ to brain wave medical systems entered into a pattern recognition computer system 

and methods and more particularly to an electroencepb- P^crably a neural network 

alograph (EEG) system and method having improved ja ^d shape is c assified as falling into one class of head 
<matial resolution 10 shapes. The electrode locations are mathematically 

2 Related Art derived in a computer system based on the known loca- 
It is presently known that the brain waves of a human lions of the corresponding electrodes on a head model 

subject, at the microvolt level, may be amplified and in the same head shape class and the size of the subject s 
analyzed using electroencephalograph (EEG) equip- h** 0 "- ctrr 
ment. Generally 19 electrodes are electrically con- 15 In another aspect of the present invention, the EEG 
nected to the scalp of the subject and the brain waves spatial resolution is improved using information about 
are amplified and displayed in analog form. The U.S. the geometry and conductivity of the skull to remove 
Pat. No. 4,736,751, incorporated by reference, describes distortion due to transmission from the brain to the 
a system using a larger number of electrodes and vari- scalp. This requires 1) Measuring the size and shape of 
ous digital computer based methods to obtain more 2° the subject's head; 2) Measuring the local skull thick- 
information from the brain wave signals. ness using an ultrasonic source and analyzing the re- 
However, when many electrodes are used, for exam- fleeted waves from the skull; 3) Estimating the local 
pie, 256 electrodes, it is time consuming to obtain their conductivity of the skull; and 4) Using these measure- 
location on the subject's scalp using the conventional ments to correct the EEG signals for distortion, 
technique of tape measurements. If the subject is a child, 25 DESCRIPTION OF THE DRAWINGS 
senile or infirm, the subject may not be able to hold still 
for the 15-30 minutes which may be required for such In the drawings 

measurements. The physical location of the electrodes FIGS. 1A-1C are diagrams of a 124-channel scalp 
is required to compare the EEG recordings, at each montage. Electrode names are based on 10-20 system 
electrode, with the information obtainable concerning 30 additional letter prefixes and numerical suffices, 
the head of the subject from other medical image sys- FIG. 1A shows the top view of the head, FIG. IB 
terns, such as MRI (Magnetic Resonance Imaging). shows the right side of the head, and FIG. 1C shows the 

The information produced from the EEG recordings j e ft side of the head, 
may be made more clear and meaningful when the brain FIG. 2 is a perspective view of the apparatus used to 
wave data is enhanced and analyzed by various meth- 35 measure the electrode coordinates. The research assis- 
ods executed in a digital computer system. However, tant touches each electrode with a magnetic position 
the presently available methods do not provide suffi- sensor which the subject is resting his head on a chin 
cient spatial resolution for some purposes. rest \^ u \\ x - m i Q a nca d support which gently restricts 

SUMMARY nead movement durin 8 the measurements. The 3-D 

40 coordinates of each electrode position are transmitted 

A functional-anatomical brain scanner with a tern- tQ the da1a co n e ction computer, 
poral resolution of less than a hundred milliseconds pj G 3 s jj 0WS tnc top views of seven representative 
measures the neural substrate of higher cognitive func- head snapeSt 

tions, as well as diagnosing seizure disorders. Electro- FI q 4 shows output of automatic edge-detection 
physiological techniques, such as electroencephalo- 45 0 f a midline sagittal Magnetic Resonance 

grams (EEGs), have the requisite temporal resolution Image obtainec i w j t h Tr-600 and TE-20 ms. The 

but their potential spatial resolution has not yet been lmes shown correspond to locations of local maxima in 
realized. Progress in increasing the spatial detail of the gradient of the image intensity. These contours pass 
scalp-recorded EEGs and in registering their functional the points at which the local average intensity 

information with anatomical models of a patient's brain 50 fe c ^ aii ^ n ^ most ra pidly. 

has been inhibited by the lack of a convenient means of nG fi show$ automated construction of finite ele- 
rapidly placing many electrodes on the scalp, determin- mem$ withm brain> skul , and ^ vo i umeS( ter- 
ing their location, determining the local thickness and rahedra] elements generated in an outer ring in brain, 
conductivity of the scalp and skull, incorporating this nG 6 fe a perspeclive view of tne flexible measure- 
information into a mathematical model, ^ and using .the 55 mem ^ ^ . piezo ^ lectric ribbons embedded 
mathematical model to enhance the spatial detail or the fa ere fa 
EEG signals. The three-dimensional positions of each 

electrode and the shape of the head are rapidly deter- DETAILED DESCRIPTION 

mined the local thickness and conductivity of the skull L ££ G Recording Method 

and scalp are determined, and this information is used to 60 

reduce blur distortion of scalp-recorded EEGs, in effect A. Nomenclature for 256 Electrode Placements 
mathematically placing the electrode just above the A nomenclature is defined for the extra equidistant 
surface of the brain. coronal rows added between the original rows of the 

In accordance with one aspect of the present inven- conventional International 10-20 System as well as the 
tion, a method is presented to improve the spatial reso- 65 extra equidistant electrodes added to fill in the spaces in 
lution of electroencephalograph (EEG) medical im- each row (FIGS. 1A-1C). For example, the preferred 
ages. The physical location of each of the electrodes in number of electrodes is 64-512 and most preferred are 
a large set of electrodes placed on the head of a subject 256 scalp electrodes. The letter M a" is prefixed to the 
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conventional 10-20 row name to indicate a position head. The patient's head shape is classified by: (1) mea- 

amerior to the existing one. For example aO is anterior suring hat stretch of stretch hat 20 haying elastic bands 

occipital, aP is anterior parietal, aC is anterior central 21 at several points over the head using piezo-electric 

and aF is anterior frontal. Two additional numbers sensors 22 as ribbons embedded within the hat (see FIG. 

appear after the location number of the 10-20 5 6; (2) scaling the vector of stretch measurements, 

positions. These two numbers (0-9) indicate the propor- termed the "stretch vector/* to have unit magnitude; 

tional distance to the next anterior and medial \0% and (3) feeding this vector to a neural network or other 

electrode position, respectively. For example, an elec- mathematical computer pattern classifier that has been 

trode halfway between P3 and aPl would be named trained to make the desired classification. The scaling 

P355, while an electrode halfway between P3 and PI 10 factor is the ratio between the stretch vector magni- 

would be named P305. With the original nineteen elec- tudes of the subject and that determined from a canoni- 

trodes of the 10-20 System, the typical distance be- cal head of average size. For each electrode montage 

tween electrodes on an average adult male head is about provided with the hat and each canonical head shape, 

6 cm; with 124 electrodes, the typical distance is about accurate measurements of electrode position are made 

2.25 cm; with 256 electrodes, the typical distance is 15 on an average sized model of that head. Appropriately 

about 1 cm. scaling these electrode positions gives a reasonably 

B. Smart EEG Hat accurate estimate of actual electrode position and has 

In most routine clinical EEGS, electrodes are pre- the advantage that no additional time is added to the 

pared and positioned on the head individually, while in EEG recording session. 

some labs commercially available electrode caps with 20 Another alternative procedure is to place a flexible 
up to 32 built-in electrodes are used. Extending this grid of electric field sources in close proximity to the 
idea, we developed a more efficient system for record- head and then measure the potential induced in each 
ing up to 124 or more EEG channels with little or no electrode when each node in the grid is activated. For 
preparation of the scalp, as described in U.S. Pat. No. example the nodes (electric field sources) may be reso- 
5,038,782, incorporated by reference herein. 25 nam circuits or switch points which form a grid on a 
C Measuring Electrode Positions flexible cap. The relative position of each node is calcu- 
Traditionally, individual electrodes are placed ac- lated by interpolation from the pattern of induced volt- 
cording to measurements taken with a tape measure. ages. The grid may be embedded in a hat which fits 
While this has been sufficient for qualitative examina- over the hat containing the recording electrodes or the 
tion of strip chart tracings with up to 21 electrodes, 30 recording electrodes and field source may be combined, 
greater precision and a method using less time are This has the advantage of increased precision with no 
needed when many more electrodes are used and when significant increase in recording time, 
one wishes to relate a recording position on the scalp to D. Software System for EEG Analysis 
the underlying brain anatomy. Although commercial EEG computer systems have 
In accordance with the present invention, the posi- 35 improved during the past few years, some important 
tion of each electrode on a subject's head is measured limitations still restrict their utility, including: (1) sys- 
with a probe 12, such as the 3-Space Isotrak digitizer by tern capacity that is inadequate for recording large num- 
Polhemus Navigation Sciences, that has coils for sens- bers of channels (124 or more) and collecting and stor- 
ing the three-dimensional position of the probe tip with ing large databases; and (2) spatial sampling and analy- 
respect to a magnetic field source attached to the head 40 ses that are inadequate, with too few channels, a lack of 
support 10 (FIG. 2). Adjustable guides 11 built into the effective means to reduce volume conduction distor- 
head support hold the subject's head comfortably in tion, a lack of cross-channel analyses (e.g., crosscovari- 
place while the measurements are made. A menu-driven ance and crosspower, correlation and coherence,), and 
program on a control digital computer 13 is used to a lack of means to investigate the relationship between 
select electrodes to be measured and display the digi- 45 neuroelectric data and cortical anatomy and physiology 
tized position of each electrode on a two-dimensional revealed by MRI and PET imaging technologies, 
projection display. Position measurement is accurate to 1. Data Collection 

better than 3 mm [RMS]. While this procedure provides One or two computers are used for stimulus prescnta- 

greatly increased precision over the use of a tape, it can tion and data collection. One computer is used to pres- 

take fifteen minutes to measure the positions of 124 50 ent stimuli and gather behavior response data from the 

electrodes during which time the patient must keep his subject while the other collects physiological data and 

head still in the head support device. Such lack of mo- controls the first computer. When a very fast computer 

tion for that time period may be difficult or impossible is used, or not many channels are required, a single 

for very young or infirm patients. computer can be used to perform all operations. A con- 

An alternative procedure is to attache the transmitter 55 trol software system runs the two computers, presents a 

for the Isotrak digitizer to the head with an adjustable variety of visual, auditory and somatosensory stimuli to 

elastic band, eliminating the need for the head support the subject according to flexible task protocols, and 

device. This alternative is more comfortable for the digitizes up to 256 channels of evoked brain wave physt- 

patient, but does not reduce the time required to mea- ological data. Up to 128 EEG traces can be monitored 

sure each electrode. 60 in real time as data are collected. Most parameters of an 

In accordance with the present invention, in another experiment can be altered via a menu-driven interactive 

alternate procedure, electrode positions are determined display. A calibration module numerically adjusts the 

by classifying a subject's head shape as one of seven gain of all channels according to the magnitude of a 

canonical head shapes (FIG. 3), determining a scale calibration signal which is injected, under computer 

factor which is the size of the subject's head relative to 65 control, into each electrode. Another module detects 

an average sized head of that canonical shape, and esti- gross artifacts and color-codes contaminated data such 

mating electrode positions as scaled values of the stan- as eye movements, gross head and body movements and 

dard electrode positions for that particular canonical bad electrode contacts. Stimulus, behavior and physio- 
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logical data are stored on hard disk according to a self- 
recording data description language, and archived to 
optical disk or magnetic tape. The data are immediately 
available to researchers at their desks using either re- 
mote terminals, or UNIX workstations via Ethernet and 
Network File System. 
2. Data Analysis 

For subsequent data analysis the system also has a 
number of other functional improvements over its pre- 


In order to visualize how the individual MRI slices fit 
together to form a three-dimensional object, a means is 
needed to juxtapose an arbitrary number of slices from 
sagittal, coronal and horizontal scanning planes into a 
composite image, and view the juxtaposed slices from 
an arbitrary viewing angle. 

C. Contour Extraction and External Surface Render- 


ing 


In order to visualize the external three dimensional 


decessors including: (1) least-squares and 3-D spline 10 surface of the brain and head, a means is needed to 


Laplacian derivation estimation and spatial deblurring 
using finite element or integral methods to reduce vol- 
ume conduction distortion; (2) digital filters with user- 
specified characteristics; (3) time series analysis includ- 
ing spectral analysis, Wigner Distributions (Morgan and 
Gevins, 1986), and event-related co variance analysis 
(Gevins et aL, 1987, 1989; Gevins and Bressler, 1988); 
(4) neural-network-driven pattern recognition to ex- 
tract optimal or near-optimal subsets of features for 


compute the curved surfaces from the information con- 
tained in the set of two-dimensional MRI slices. Scalp 
and cortical surface contours are extracted from MRIs 
using two image analysis methods we developed. The 
15 first technique used intensity thresholding, which in- 
volves extracting contours along which the image in- 
tensity is equal to a defined threshold. This technique is 
useful for extracting the scalp surface contour which 
can be discerned easily from the black (approximately 


recognizing different experimental or clinical categories 20 zero intensity) background of the image. The technique 


(Gevins, 1980, 1987; Gevins and Morgan, 1988); and (5) 
anatomical modeling to construct 3-D mathematical 
models of the brain and head from MRI to CT brain 
scans. Four on-line, interactive subsystems are used to 
examine and edit data for residual artifacts (on an indi- 
vidual channel basis if desired), sort data according to 
stimulus, response or other categories, perform explor-. 
atory data analysis, and produce three-dimensional 
graphics representations of the brain and head. 

II. Brain Modeling Methods 

Since commercial magnetic resonance image (MRI) 
analysis packages lack essential features for functional- 
anatomical integration studies, there is described below 


25 


30 


has six steps: 1) for each pair of adjacent contours, find 
the point on the second contour closest to the first; 2) 
calculate the two distances from each of these points to 
the next point on the adjacent contour; 3) make a trian- 
gle using the point with the shortest distance, advancing 
on that contour; 4) repeat steps two and three until all 
points are exhausted; 5) repeat steps two to four in the 
reverse direction along the contours; and 6) piece to- 
gether the "best" results. 

The second contour extraction method involves dif- 
ferential intensity analysis. Using this technique, con- 
tours that separate image regions with different local 
average intensities are extracted. Resulting contours 
pass through the points in the image at which the local 


algorithms and a soft ward outline to produce 3-D brain 35 average image intensity is changing most rapidly. This 


models suitable for functional localization studies. Visu- 
alization softward permits construction of 3-D compos- 
ites of multiple 2-D image planes, as well as 3-D surface 
rendering based on surface contours. Since generating 
surface contours manually is laborious and subject to 40 
error, we have largely automated the procedure. We 
have also automated the alignment of the digitized EEG 
electrode positions with the scalp surface contours, 
which is a critical first step in a functional-anatomical 
analysis. The MRI contour information is used to pro- 45 
duce mathematical finite element or surface integral 
models of the brain and head suitable for equivalent 
dipole source localization and scalp EEG deblurring 
procedures. 

A. EEG-MRI Alignment Procedure 50 
In order to visualize the brain areas underlying EEG 
electrodes, a procedure is needed for aligning scalp 
electrode positions with the MR Images. In our proce- 
dure, x, y, z translation and x, y, z axis rotations are 
computed iterative! y to align the digitized positions of 55 
the EEG electrodes with the MRI data. This is done for 
each electrode by finding the distance to the closest 
point on the scalp surface MRI contours and minimiz- 
ing the mean distance for all electrodes. With MR im 


technique requires no a priori intensity threshold value, 
and is useful for extracting the cortical surface contour 
which has a less well-defined image intensity value 
throughout an image than the scalp surface contour. 
The first -order and second -order partial derivatives of 
the image are estimated using 2-D filters, and these 
derivatives are used to locate the local maxima in the 
gradient of the image intensity. Highly computationally 
efficient filtering techniques have been developed for 
the estimation of the partial derivatives of the image 
(Algazi et al., 1989). FIG. 6 shows an example of con- 
tours corresponding to locations of local maxima in the 
gradient of the image intensity, which trace sulci. 

D. 3-D Cortical Internal Surface Image Model 3-D 
surface models of the external convexity of the cerebral 
cortex are not sufficient to visualize and computation- 
ally model the cortical surface within fissures and sulci. 
As a result, we developed an algorithm to model cubic 
volume elements (voxels). The faces of the voxels lie in 
the horizontal, coronal, and sagittal planes for which 
MRI data have been obtained. The MRI data are 
mapped onto the faces of the voxels to obtain a 3-dimen- 
sional display of the image data. The image planes used 
are averaged horizontal planes that lie halfway between 


ages that have a 3 mm inter-slice spacing, the mean 60 the horizontal planes in which MRI data have been 


error distance is usually better than 2 mm. This is more 
accurate and less subjective than alignment procedures 
that use skull landmarks such as the nasion and pre- 
auricular points located visually in the MR images. 
FIG. 1A shows the electrodes displayed schematically 65 
on a scalp surface reconstructed from horizontal MR 
scalp contours as described below. 
B. 3-D Composite MRI Displays 


acquired, and coronal and sagittal image planes (that 
can optionally be synthesized from the acquired hori- 
zontal images). The spacing between image planes is 3 
mm, which yields voxels with dimensions of 3 mm by 3 
mm by 3 mm. The images have pixel dimensions of 
approximately 1 mm by I mm. 

The initial set of voxels is the set bounded by those 
voxels that lie just inside the cortical surface. To view 
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the MRI data at a slightly deeper level, a mathematical 
morphology erosion operation is used to remove the 
boundary layer of voxels, thereby exposing the faces of 
the vosels that lie one layer deeper. By eroding the 
model iteratively, anatomical structures can be tracked 5 
and 3-dimensional models of the structures can be made 
(see FIG. 7). 

. E. 3-D Electrical Brain Model: Finite Element 
Method 

A mathematical 3-D representation of the brain and io 
head are needed to most effectively deblur the EEG 
recorded at the scalp and to compute equivalent current 
dipole sources from the scalp EEG. We developed an 
efficient finite element method for this purpose. 

Maxwell's Equation, IS 


V<crVy)-?-V.J=0 in ft 


(1) 


is frequently used to study the electromagnetic field 
generated by populations of neurons. Here <r(>0) is the 2Q 
conductivity tensor at a point =(x,y,z) in fl(c.g M a 
human head), u is the electric potential at R, J is the 
electric current density at R and V is the grad vector 
operator. When fl represents a human head, the follow- 
ing boundary condition for equation (1) is obtained 25 


a-Vu n =0 on £(headsurface) 


(2) 


since electric current does not flow out of the head in 
the direction normal to the surfaces. The limited scalp 
potential measurement specifies another boundary con- 
dition. 


30 


u= U(x,y.z) onS]. 


(3) 


Where Si is a subsurface of S corresponding to the area 35 
covered by the electrodes. 

The Finite Element Method (FEM) we developed 
finds the numerical solution of equation (1) with the 
boundary condition as specified in equation (2) when 
the electromagnetic field activities of a head are mod- 40 
eled. Our implementation has the following special 
features: 

a) A procedure to generate a list of wedges within a 
head such that each wedge has an unique tissue type 
and the interfaces between two neighboring wedges 45 
are perfectly matched to guarantee a correct numeri- 
cal solution to Eq. (1). 

b) A proper numbering scheme on a chosen set of nodes 
and the use of local supporting functions to produce 

a diagonally condensed block sparse coefficient ma* 50 
trix. 

c) The use of elementwise local shape functions to rep- 
resent the local supporting functions. 

d) The use of a compressed data structure to obviate the 
need to allocate a very large memory space so that 55 
large size problems can be solved on desktop com- 
puter workstations. 

e) Use of a C-language sparse matrix solver incorporat- 
ing a temporary working array in the process of Cho- 
lesky Decomposition to speed up the process of com- 60 
puting the numerical solution. 

Advantages of this approach are: 1) it handles the 
Dirac Delta Function V J smoothly by transforming 
equation (1) into a variational form when J represents a 
dipole-like kind of source; 2) it allows modeling of the 65 
complicated geometries of different tissues within the 
head by generating finite elements using contours on 
pairs of adjacent MR images; and 3) it produces a sparse 


matrix in which many entries are zeros and thus is con- 
venient and fast to compute on a desktop computer 
workstation. For example, a brain model with 10,087 
FEM nodes can be solved in about 100 minutes on a 
SUN Sparc-1 workstation (12 MIPS) with 8 MB of 
memory (FIG. 5). 

111. Spatial Enhancement Techniques to Reduce Blur 
Distortion 

Electrical currents generated by sources in the brain 
are volume conducted through brain, cerebrospinal 
fluid, skull and scalp to the recording electrodes. The 
principal source of distortion in this process is produced 
by the skull, which has low conductance relative to the 
other tissues. Because of this spatial low-pass distortion, 
the potential distribution of a localized cortical source is 
spread over a considerable area of scalp and appears 
blurred or out of focus. For example, using a 4-shell 
spherical head model, we estimated the "point spread" 
for a radial dipole in the cortex to be about 2.5 cm. 
There are several methods to reduce this distortion 
which involve a tradeoff between the amount of infor- 
mation about head shape and tissue properties required 
versus the accuracy of the distortion correction. 

A. Laplacian Derivation 

The simplest method is the spatial Laplacian operator 
or Laplacian Derivation (LD). It involves computing 
the second derivative in space of the potential field at 
each electrode. This converts the potential into a quan- 
tity proportional to the current entering and exiting the 
scalp at each electrode site, and eliminates the effect of 
the reference electrode used during recording. An ap- 
proximation to the Laplacian Derivation, introduced by 
Hjorth (1975, 1980) assumes that electrodes are equidis- 
tant and at right angles to each other. Although this 
approximation is fairly good for some electrode posi- 
tions such as midline central (Cz). it is less accurate for 
others such as midtemporal (T5). We have developed 
two methods which produce more accurate estimates of 
the Laplacian using the actual measured electrode posi- 
tions. 

1. Planar Projection and Least Squares Solution 
The first method is based on projecting the measured 
electrode positions onto a two-dimensional surface and 
estimating the surface Laplacian with a I east -squares 
procedure. 

The objective of the LD calculation is an accurate 
estimate of 


(4) 


Assume we wish to compute the LD at position q 0 
using N surround electrodes in an arbitrary configura- 
tion. Let the voltage at position q/be designated by V,. 
Put the origin of a cartesian coordinate system at q 0 
with the z axis pointing perpendicular to the plane near- 
est the electrodes. Let (qixjiy) represent the x and y 
components of the electrode position within this plane. 

Then E(x,y)= -VV(x,y) is the electric field in the 
plane or equivalently: 


E-ds 


(5) 


01/04/2004, EAST version: 1.4.1 


5,119,816 


Where C is any curve in the plane connecting q 0 to q,. 
If we let Ex and E v designate the x and y components of 
E then the two-dimensional laplacian of the voltage is 
given by: 


(6) 


The field E can be expanded in a two-dimensional 
Taylor series (Rudin, 1976; Apostol, 1969) and ex- 
pressed in a single matrix equation: 


10 


QF~V 


(7) 


15 


In this equation the Q matrix is composed of the x and 
y coordinates electrode positions and remains fixed 
during the entire recording. The voltage vector V is 
measured at each time point and vector F composed of 
the electric field components and their derivatives is 20 
unknown. 

Multiplying the above by the transpose of Q gives 


10 

(£o f T)o) is a point in Jij-plane. For convenience, we de- 
note ax/a £ by fs6, etc. By the definition of the Laplacian 
here, f, g and h have to be the transformation functions 
for points on the £T)-plane to points on the surface S 
defined in the current xyz-Cartesian coordinates. Gen- 
eral forms of f, g and h for all points on £17- plane are 
difficult to find. However for every point (xo, yo, zo) on 
S t explicit forms of f, g and h can be found for transla- 
tion and rotation transformations of points from £17- 
plane to the tangent plane of S at (xo, yo» zo) in xyz- 
Cartesian coordinates. The explicit form of the potential 
distribution function U(x,y,z) on a scalp surface S is not 
available in EEG recordings. However a set of scalp- 
recorded potential V/(i= 1,2, • - . n) is available on a set 
of recording channels on S, (x/, yj, z/), i=l,2, . . . , n. 
Therefore the potential distribution can be computed by 
interpolation from the scalp-recorded potentials. First a 
set of basis functions <f>,<x,y,z) must be computed such 
that 


Q T QF=Q T v 


(8) 


If Q premultiplied by its transpose Q r is nonsingular 
then we can form the least squares estimate of F as 


25 


F=[Q T Q)~ l Q T V 


(9) 


If we let C be the sum of the third and fifth rows of 
[Q 7 "Q]"" , Q r then from the above equation we can esti- 
mate the laplacian LD as: 


30 


v 2 v 


40 


This method accounts for both variations in distance 
and variations in angle between electrodes in a compu- 
tationally efficient manner, since once the coefficients 
are computed based on the positions of the electrodes, 
application of the method only requires a single matrix 
multiplication for each point of a time series. Although 
this procedure produces a dramatic improvement in 45 
topographic detail, some problems remain because of 
the projection onto a two-dimensional plane, and the 
assumption that the current gradient is uniform over the 
region encompassed by the surrounding electrodes used 
to estimate the Laplacian of an electrode. 

2. Three Dimensional Spline Method 

A more accurate method, of estimating the Laplacian 
Derivation does not rely on the two-dimensional pro- 
jection. 

We compute the Laplacian of a distribution function 
U(x,y t z) (defined on a convex surface S) at its tangent 
plane attached to S at point (xo,yo»zo) where the surface 
S is defined by: 


50 


55 


and 


(12) 


(13) 


Since EEGs can be described by Maxwell's Equa- 
tion, the open solution of that equation can be used for 
the basis functions. Accordingly, we define: 


(10) 35 


((* - x/o) 2 -»- (v - yiolr + <r - j*)) 2 )* 


(14) 


We want to compute: 


such that 


an 2 


at io. "TO. 


(11) 


Here (x*), y/o. Z/o) is the center of the best fitting 
sphere of radius 1 to the point (x,-, y; f z,) on the surface 
S. Algorithmically we choose this center point 1 unit 
from the i-th channel position (x„ y,\ z,) in the opposite 
direction of the normal orientation of surface S at (x/, y,\ 
z/). The basis function <J>, is a local support function and 
is smooth enough in the neighborhood of (x/, yy, z,*). For 
a given channel list (x/, y„ z/) and its corresponding 
scalp potential data set V/, (i=l,2 n). we 

1) construct a polygon montage spanned by these 
channels, that defines the surface S; 

2) estimate the normal directions for every channel 
position by averaging the normal orientations of 
polygons surrounding this channel position; 

3) construct a set of basis functions <f>/(x,y,z) and solve 
the weights a/ through systems of equation in (13) 
by the Gaussian Elimination Method; 

4) for every point on S, the desired surface tangent 
laplacian can be obtained by equations (11) and 
(12)- 

Equation (12) suggests that the desired surface tan- 
gent laplacian of U(x f y,z) can be computed by taking 
the surface tangent Laplacians of <J>/(x,y,z) for every 
60 point (x,y,z) on S, and then linearly combining them 
with weights a/. 
3. Finite Element Method 

A better improvement in distortion correction is pos- 
sible by using a Finite Element Method to represent the 
true geometry of cortex, cerebrospinal fluid, skull and 
scalp, and estimates of the conductivities of these tis- 
sues, to model the potential activity described by Max- 
well's equations in the scalp and skull layers. This aJ- 
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lows estimation of the potentials which would actually from u/* and U{x,y.z) which correspond to the pre- 

be recorded on the surface of the brain. This can be dieted and measured potentials, respectively, 

done without introducing an arbitrary model of the Equation (19) puts forward an alternative to the 

actual number and location of sources because, regard- above approach. For each ui, we can compute the corti- 

less of where they are generated in the brain, potentials 5 cal potential U3. If there are a large number of recorded 

recorded at the scalp must arise from volume conduc- potential observations, such as hundreds of sampled 

tion from the cortica] surface through the skull and time points, the following process will reduce the calcu- 

scalp lation time dramatically. Suppose e,- is a vector of the 

We use Maxwell's equation and our Finite Element same dimension as U] with all its components as zero 

Method to relate the scalp-recorded potentials with the 10 except 1 at i-th position. Then 
cortical potentials without introducing an explicit 

source function. For this application, the region of in- ni 

terest is limited to the scalp volume and skull volume. u 1 = t*i. ^ • • • • = * w 
The difference between the inner surface of the skull 

and the outer convexity of the conical surface is as- 15 lf we soIve (he uain equat j on (19 ) w j t h Ul replaced by 

sumed to be small enough to be neglected. We also e / and denote the resulting solution vector by p f( then for 

assume that the potential activity in the region of inter- an ^ bilT2kTy U /=(si,s 2 , . - . ,s fll )', the corresponding solu- 

est is described by Maxwell's equations, that the bound- t j on j s; 
ary condition on the air-scalp surface is stated by equa- 
tion (2), i.e., current cannot flow out of the scalp, and 

that the boundary condition on the cortical surface is U) w J sp . 

stated by the following equation which describes the I 
cortical surface potentials: 

which is a simple linear combination of ni vectors and 

w=C{x.y,z) on s 2 . (15) ^ 5 therefore much faster to compute. 

4. Surface Integral Method 

We also assume that there is no generating source An alternative to performing deblurring (or dipole 

within the scalp or skull which means the function J in source localization) through using numerical solutions 

equation (1) is zero. Although of course there is no to the differential form of Maxwell's equations is to use 

general solution to the "inverse problem" because of the 30 numerical solutions to the integral form of Maxwell's 

assumed boundary conditions there is a unique solution equations. 

of equation (1), which may be obtained as follows. Ap- The following equation is used for the source local- 
plying the Finite Element Method to equation (1) with ization problem: 

boundary conditions in (2) and (15), the following ma- ^ (20) 


trix vector relation is obtained: 


a — 

■ R Rr- 


IS. 


Au-f (16) <™«> u - jZt) s .^ v '' ou ' Rl */ 

Here u is a vector of dimension n which is a numeri- where denotes conductivit iois a inl on the outer 
cal approximation to the analytical potential distnbu- surface ^ ^ u i$ ja , g h ^ vcc(or from 
tion function u m equation (1). The value of n corre- & source to ; £ ^ kude / 3 is the moment of 
sponds to the total number of vertices on all the finite fl , e m J £ thc numbcr of imcrfaccs> s . ^ 
elements in the region of interest. Since we are assuming tween ion$ ^ conductivity . n . * the unit 

that no generating sources are present in the skull and 45 yeclor frQm thfi ion t0 the QUt ion; and S is thc 
scalp J=0 and therefore f=0. We then decompose u vector from xo to the point of integration with R denot- 
mto three sets, the potentials which correspond to the . ^ unh vectQr and R denoling ^ magnitude. A 
electrodes on the scalp, the potentials which corre- numer ical solution of this equation involves obtaining a 
spond to the cortical surface, and the potentials m the discrete approximation to the ime g ra J. Boundary ele- 
rest of the region, and denote them by ui, uj and u 2 5Q mcnt mcthodscanbc do this. An alternative is to 

respectively, If we decompose the matrix A corre- w efficient integration methods that have been devel- 
spondmgly, then equation (16) becomes: opcd for thc sphcrc md methods of computational ge- 

ometry for projecting a non-spherical surface on a 
0?) sphere. 

55 One advantage of using the integral method for 
source localization is that no nodes are close to the 
sources, a locus where potential fields vary a great deal. 
This is a problem for the finite element method. 
Solving ui with respect to u 3) we get ™e following equation is used for the deblurring 

60 problem: 

W\\-A\iAn-*A2x)v^{AxiAn- X Al*-A\iiu> (18) 

and solving U3 with respect to ui, we get 2 <ow - ov.w) — - — u = 2*<r76u 

'° * J Si - XQ K ~ 

where the notation is that for Equation (20). After regu- 
For a given Gjt, U3 fc is known. Then ui* can be com- larization, the numerical integration techniques used for 
puted via equation (18) and the residual is calculated source localization are applied. An advantage of this 
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formulation when applied to deblurring is that it ac- 
counts for the geometry of the whole conductor 
whereas the finite element technique only accounts for 
the portion of the conductor that includes the elec- 
trodes on the scalp and the tissue immediately below. 
5. Estimating Scalp and Skull Thickness for Models 
In the finite element or surface integral methods, the 
local geometry of scalp and skull are each explicitly 
represented. If an MRI or CT brain scan is not avail- 
able, an alternative method of measuring tissue thick- 
ness can be used based on low-frequency ultrasound 
technology. 

For a low-resolution representation of skull thick- 
ness, the "end echo" in standard echoencephalography 
is used to estimate the thickness of the skull. This re- 
quires placing the ultrasonic probe at a position on the 
opposite side of the head from the place of desired mea- 
surement, normal to the scalp at the probe, and normal 
to the skull at the place of skull thickness measurement 
on the opposite side of the head from the probe. Skull 2 o 
thickness at each of the major cranial bones is measured 
in the following locations: 


10 


15 


crania) bone 


position probe position 


left temporal T3 near T4. directly above the ear 

right temporal T4 near T3. directly above the ear 

frontal (left) aFl-aF3 roughly 1 cm below 02 

frontal (right) aF2-aF4 roughly I cm below Ol 

left occipital Ol .2-.5 cm righi of FPZ 

right occipital 02 .2-. 5 cm left of FPZ 

left parietal aP3 I cm below M2. with upward till 

right parietal aP4 1 era below Ml ■ with upward tilt 


25 


30 


The end -echo consists primarily of the echo due to 
reflection from soft-tissue to inner skull interface and 
the echo due to outer skull to scalp interface (followed 35 
by a very small scalp-to-air echo). The delay prior to 
the end-echo varies, according to scalp location and 
inter-subject variation, roughly from 90 to 120 micro- 
seconds. The delay between the echos from either sur- 
face of the skull may vary from 1.7 to 5 microseconds. 40 
The second echo is attenuated with respect to the first, 
and requires a gain increase of 2.5 dB per MHz (typical) 
to be comparable in size to the first. Pulse widths of 0.25 
to 1 microseconds are used with burst frequencies of 2 
to 10 MHz, adjusted to produce the best discrimination 45 
between the two echos. 

As necessary, passive receptors are used to scan the 
area of the intended skull thickness measurement to 
determine the accuracy of aim of the transmitted beam. 
When necessary, an average is taken over a small clus- 
ter of locations (3 to 6) near the listed probe site, in 
order to average over the variations in thickness for a 
given bone near the target area. In those instances when 
the aimed beam cannot be focused well on the parietal 
bones, the average thickness of these bones is estimated 
from normalized data on the relationship between pari- 
etal bone thickness and the average thickness of the 
other bones. 

For finer spatial resolution in measurement of skull 
thickness, a near side detector can be used. It has two to 
four transmitters adjacent to the transmitter/receiver 
probe, spaced far enough away so that the critical angle 
(27 degrees) prevents energy from the adjacent trans- 
mitters from entering the skull at an angle that could be 
reflected back to the probe. 

The phase and time delay between the adjacent trans- 
mitters and the probe is adjusted to minimize the rever- 
beration at the probe. Thus, at the probe site the adja- 


cent transmitters are contributing energy solely through 
reflections from the outer skull surface, and in a way to 
minimize reverberation at the probe. Only energy trans- 
mitted by the probe itself can enter the skull and be 
reflected back to the probe. 

The thickness of the near skull can then be measured 
either by subtracting independent measures of the re- 
flection time from the outer skull surface (adjacent 
transmitters inactive) and reflection time of the echo 
from the inner skull surface (reverberation minimization 
active), or by interferometry. The interferometry op- 
tion sweeps the sound frequency (readjusting reverber- 
ation minimization for each frequency), and determines 
the change in frequency required for one oscillation in 
the variation in reflected energy due to interference of 
the inner and outer reflections from the skull. The mea- 
sured thickness of the skull is one half the speed of 
sound in the skull divided by the change in frequency 
required to move through one "interference fringe". 
6. Estimating Conductivities 

Up to now, estimates of conductivity for the scalp, 
skull, CSF and brain have been from published sources 
(Geddes & Baker, 1967; Hosek, 1970). More accurate 
values are needed to realize the full benefit of the FEM 
head model. This can be done by solving for the ratio of 
skull-to-scalp and skull-to-brain conductivities using the 
scalp potential distribution produced by a single com- 
pact population of neurons, represented as an equivalent 
current dipole, whose location is approximately known. 
In one procedure, a steepest-descent, non-linear, least- 
squares fit between the measured potentials and the 
forward solution of an equivalent current dipole in a 
three-sphere model is performed (Fender, 1987). The 
main unknown parameters are the ratios of conductivi- 
ties of the tissues, as well as the dipole strength (which 
is not of interest here). The availability of reasonable 
initial guesses, as well as upper and lower bounds, from 
existing experimental data result in realistic solutions for 
the unknown parameters. A highly localized source 
suitable for this purpose is the somatosensory evoked 
potential peak at 20 msec produced by mild stimulation 
of the contralaterial median nerve at the wrist. This 
potential, negative posterior to the central sulcus and 
positive anterior to it, is generated by a single compact 
population of neurons oriented tangentially in the hand 
area of somatosensory area Zb and located about 2.5-3.0 
cm below the surface of the scalp. Published values of 
the amplitude of the cortical potential and the ratio of 
cortical to scalp potentials are also available (Dinner et 
al., 1983). The global conductivity parameter is then 
adjusted for each of the finite elements in the skull using 
local skull thickness measures determined from the low 
frequency ultrasound, MRI or other means. 

An alternative method for estimating tissue conduc- 
tivity is based on measuring the reverse electromotive 
force produced by eddy currents induced in the tissue 
by an applied external field generated by a coil placed 
60 against the head of each patient. The magnitude of this 
emf depends on the conductivity of the tissue and the 
frequency of the magnetic field, i.e. higher conductivity 
means more eddy current is induced which gives a 
larger reverse emf. Very low frequencies induce eddy 
currents in all tissues whereas high frequencies induce 
currents only in the tissues closest to the coil. Scanning 
a spectrum of frequencies yields data from which the 
conductivities of each tissue type can be computed. 
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We have also used FEM to: 1) localize single equiva- 
lent dipolc sources of recorded scalp potentials; 2) ob- 
tain the potential distribution within the whole head 5 
given the localized source. We choose a specific func- 
tion J and solve the corresponding potential distribution 
u using FEM such that the difference between u on Si 
and U(x,y,z) defined on Sj is minimized. The following 
source localization procedure (SLP) outlines an itera- 10 
tive process used to find J and u. For a chosen Jo and 
K=0, 1, 2 repeat the following: 1) find the numeri- 
cal solution ujt of equation (1) on CI coupled with the 
boundary condition in equation (2) using FEM; 2) mea- 
sure the difference between u* and U(x,y,z) on Si and 15 
choose Jjt+i; 3) compute the distance between Jjl-J* + |. 
If the distance is less than some tolerance, stop the loop. 
We applied the SLP to the case that J represents a single 
dipole. We used the Simplex algorithm (Nelder and 
Mead, 1965) for computing J*+i. 20 
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I claim: 

1. The method of measuring the positions of a plural- 
ity of electroencephalographic (EEG) electrodes re- 
movably electrically connected to the head of a subject 
by: 

(a) measuring the head of the subject by measuring a 
plurality of standard cross-directional lines on the 
head to obtain a set of head measurements; 

(b) entering the set of head measurements into a com- 
puter system which automatically determines the 
classification of the head shape according to a plu- 
rality of predetermined head shape classes; 

(c) comparing the subject head size to a predeter- 
mined head model of the same class in which the 
position of electrodes corresponding in position to 
the electrodes to be measured have been deter- 
mined; and 

(d) determining the position of the electrodes con- 
nected to the head of the subject by scaling the 
measurements of the electrode positions on the 
head of the subject in accordance with the deter- 
mined electrode positions of said head model. 

2. The method of claim 1 wherein the measurements 
of step (a) are obtained by using a stretchable hat having 
piezoelectric ribbons embedded therein. 

3. The method of claim 2 wherein the measurements 
are stretch measurements obtained from stretching of 
the ribbons and including the steps of obtaining a vector 
of said stretch measurements and scaling said stretch 
measurements. 1 

4. The method of claim 1 wherein said computer 
system is a pattern recognition neural network com- 
puter. 
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